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1. Introduction 

The concept of scale-invariance is central to the modern understanding of critical 
phenomena in and out of equihbrium. Its exploitation through the renormalization 
group has in particular led to the recognition of universal critical exponents and scaling 
functions which describe the behaviour of physical observables, see e.g. P and references 
therein. Here we are concerned with the slow dynamics of systems brought rapidly to 
their critical point and/or into a phase with more than one thermodynamically stable 
state. Such a kind of behaviour is typical for glassy systems but also occurs in simple 
magnets with a purely relaxational dynamics which were quenched from a disordered 
state to a final temperature T < Tc, where Tc > is the critical temperature. For the 
latter, it is now understood that the dynamics is governed by a single time-dependent 
length-scale L = L{t) ~ t^^^ for t sufficiently large and where z is the dynamical 
exponent j2j. As an example, consider simple magnets relaxing towards equilibrium. 
For phase-ordering kinetics (T < T^), the Bray-Rutenberg theory shows that dynamical 
scaling together with the assumption of a Porod law for the time-dependent structure 
factor predicts the value of 2; whereas for T = the value of z is computed from 
critical (equilibrium) dynamics. More recently, it has been understood that the study 
to two-time observables provides further and deeper insight, in particular the ageing 
behaviour is made explicit through the breaking of time-translation invariance. The 
challenge is now to find the values of the associated non-equilibrium (ageing) exponents 
and also the form of the scaling functions, see below for the precise definitions. 

A common way to study this problem is through a Langevin equation which should 
describe the dynamics of a coarse-grained order-parameter. This may be turned into 
a field-theory and renormalization-group methods then allow to extract values of these 
exponents, in quite good but not perfect agreement with the results of direct numerical 
simulations ^ . On the other hand, the resulting predictions for the scaling functions 
appear to be far from the numerical results, see Hj . An alternative approach seeks 
to extend dynamical scaling to a larger group of 'local' scale-transformations jH], see 
[H] for a recent review. In the framework of a de Dominicis-Janssen type theory 
[Tin ITT] , the effective action S = S[(j),(f)] is given in term of the order-parameter 
field and its associated response field 0. Furthermore, for systems in contact with 
a thermal bath such that detailed balance holds one always has the decomposition 
S[(f),(j)] = S'o[0, 0] + Sb[(j)] into a 'deterministic' part Sq which can be derived from the 
Langevin equation when all noise terms are dropped and the 'noise' term 5*6 [0] which 
depends only on the response function ^2] • The form of response functions can then be 
found from the requirement of covariance under the group of local scale-transformations. 
As we shall explain below, correlation functions can be reduced to certain integrals of 
higher, multipoint response functions ^21- This approach yields the form of the scaling 
functions whereas the exponents are treated as parameters whose values have to be 
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suppliedf and reproduce perfectly the results of both analytical and simulational studies 
of many common spin systems undergoing phase-odering kinetics where z = 2.% 

It is an established fact that the basic Langevin equation for the order-parameter 
does not admit any symmetries beyond dilatations and (space-)translations, see ^1] 
for a recent discussion. However, it has been shown that at least for simple magnets 
it is enough to concentrate on the dynamical symmetries of the deterministic part of 
the Langevin equation only, as given by the action 5*0 . In particular, Schrodinger- 
invariance of that deterministic part is sufficient to be able to derive the two-time 
correlations C(t, s) and two-time response functions R(t,s) explicitly J^. There is 
an exact agreement for systems such as the spherical model, the XY model in spin-wave 
approximation or the voter model which are all described by a linear Langevin equation. 
Good agreement with simulations of Ising, Potts and XY models models was found as 

well nmniiin!. 

In this paper, we extend the treatment of local scale-invariance to ageing systems 
with a dynamical exponent z = 2 but without detailed balance. Working with a de 
Dominicis-Janssen type theory, we find again a decomposition S'[0, 0] = So[(f), (j)]+Sb[(p, 4>] 
into a 'deterministic', Schrodinger-invariant term 5*0 and 'noise' terms, each of which 
contains at least one response field more than order-parameter fields (explicit expressions 
will be given in sections 2 and 3). Then the Bargman superselection rules which follow 
from the Galilei-invariance of 5*0 are enough to establish that again the two-time response 
function is noise-independent and the two-time correlation function can be reduced to 
a finite sum of response functions the form of whom is strongly constrained again by 
the requirement of their Schrodinger-covariance. These developments provide further 
evidence for a hidden non-trivial local scale-invariance in ageing systems which manifests 
itself directly in the 'deterministic' part (see |jl7| for the construction of Schrodinger- 
invariant semi-linear kinetic equations) but which strongly constrains the full noisy 
correlations and reponses. 

We test the present framework of local scale-invariance in two exactly solvable 
systems with a non-linear coarse-grained Langevin equation. A convenient set of models 
with a non-trivial ageing behaviour is furnished by the bosonic contact [18^ and pair- 
contact processes both at criticality. These systems are defined as follows. Consider 
a set of particles of a single species A which move on the sites of a hypercubic lattice in d 
dimensions. On any site one may have an arbitrary (non- negative) number of particles. § 

f This is close in spirit to the treatment of equilibrium phase transitions through conformal invariance, 
which fixes the form of the n-point correlators in terms of the scaling dimensions of the scaling fields 
0. Furthermore, those exponents can be determined exactly in 2D from symmetry considerations 
(i.e. representation theory of the Virasoro algebra) alone since the conformal symmetry is infinite- 
dimensional in that case. 

I For non-equilibrium critical dynamics where 2 7^ 2 in general one also has a good match with 
numerical data for the response functions in direct space but systematic differences may appear in 
momentum-space calculations, e.g. in the 2D Ising model for t/s < 10 

§ This property distinguishes the models at hand from the conventional ('fermionic') contact and pair- 
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Single particles may hop to a nearest-neighbour site with unit rate and in addition, the 
following single-site creation and annihilation processes are admitted 

mA — ^ (m + 1)A , pA — ^ {p — i)A ; with rates /i and A (1) 

where £ is a positive integer such that \i\ < p. We are interested in the following special 
cases: 

(i) critical bosonic contact process: p = m = 1. Here only £ = 1 is possible. 
Furthermore the creation and annihilation rates are set equal fi = X. 

(ii) critical bosonic pair-contact process: p = m = 2. We fix £ = 2, set 2X = fi and 
define the control parameter || 

(2) 

The dynamics is described in terms of a master equation which may be written in a 
hamiltonian form dt\P{t)) = —H\P{t)) where \P{t)) is the time-dependent state vector 
and the hamiltonian H can be expressed in terms of creation and annihilation operators 
a(a:;,t)^ and a{x,t) jSOlElESl- It is well-known that these models are critical in the 
sense that their relaxation towards the steady-state is algebraically slow [IHl dl 121] ■ In 
particular, the local particle-density is p{x,t) := {a{x,t)). Its spatial average remains 
constant in time 

Jdxp(x,t) = J dx {a{x,t)) = po (3) 

where po is the initial mean particle-density. We are interested in the two-time connected 
correlation function 

G{r; t, s) := {a{x, t)a{x + r, s)) — pi (4) 

and take an uncorrelated initial state, hence G{r] 0, 0) = 0. The linear two-time response 
function is found by adding a particle-creation term Ylx^i^'^) {(^K^^ 

t) - 1) to the 

quantum hamiltonian H and taking the functional derivative 



5h{x, s) 



(5) 

h=0 



We have previously analyzed these quantities in the scaling limit where both t, s as well 
as t — s become large with respect to some microscopic reference time. The results are 
as follows [23- consider the autocorrelation and autoresponse functions, which satisfy 
the scaling forms 

G{t,s):=G{0;t,s) = s-'fGit/s) (6) 
Rit,s):=RiO;t,s) = s-'--fnit/s) (7) 
contact processes whose critical behaviour is completely different. 

II If instead we would treat a coagulation process 2 A — > A, where i = I, the results presented in the 
text are recovered by setting X — fi and a — n/D. 
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bosonic contact process 


bosonic pair-contact process 


a < Oc 


a = 


a 


2 ^ 


2 ^ 


- 1 

2 


b 


2 


d _ 1 
2 ^ 


if 2 < ci < 4 
f-2 if ci>4 



Table 1. Ageing exponents of the critical bosonic contact and pair-contact processes in 
the different regimes. The results for the bosonic contact process hold for an arbitrary 
dimension d, but for the bosonic pair-contact process they only apply if d > 2, since 
ac = for d < 2. 







fciy) 


contact process 






pair 
contact 
process 


a < ac 


d>2 






a = ac 


2 < < 4 






d > 4 




+ l)-f+2 _ _ l)-f+2 + _ 4)(^ _ l)-f+l 



Table 2. Scaling functions (up to normalization) of the autoresponse and 
autocorrelation of the critical bosonic contact and bosonic pair-contact processes. 



where the values of the exponents a and b are hsted in table ^ Here the critical value 
Uc for the pair-contact process is explicitly given by [T^ 

— = 2 / du (e-^"Jo(4M))'' (8) 
Q^c Jo 

where /q is a modified Bessel function. The dynamical behaviour of the contact process 
is independent of a. For the critical bosonic pair-contact process, there is a clustering 
transition between a spatially homogeneous state for a < ac and a highly inhomogeneous 
one for a > Uc where dynamical scaling does not hold. These two transitions are 
separated by a multicritical point at a = ac- Since our models do not satisfy detailed 
balance, there is no reason why the exponents a and b should coincide and our result 
a ^b for the bosonic pair-contact process is perfectly natural. 

While the scaling function = {y — has a very simple form, the 

autocorrelator scaling function has an integral representation 

fciy) = Go f'de e^-\y + 1 - 2dY''^ (9) 

where the values for a and b are given in table Q and Qq is a known normalization 
constant. The explicit scaling functions are listed up to normalization in table El |24j . 
In this paper, we shall study to what extent their form can be understood from local 
scale- invariance. 
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This paper is organized as follows. In section 2 we treat the bosonic contact process 
in its field-theoretical formulation. The action is split into a Schrodinger-invariant term 
5*0 and a noise term Sb and we show how the response and correlation functions can 
be exactly reduced to certain noiseless three- and four-point response functions. In 
this reduction the Bargman superselection rules which follow from the Schrodinger- 
invariance of 5*0 play a central role. These tools allow us to predict the reponse- 
and correlation functions which will be compared to the exact results of table |21 In 
section 3 the same programme is carried out for the bosonic pair-contact process but 
as we shall see, the Schrodinger-invariant term 5*0 of its action is now related to a 
non-linear Schrodinger equation. The treatment of this requires an extension of the 
usual representation of the Schrodinger Lie-algebra which now includes a dimensionful 
coupling constant. The construction is carried out in appendix A. The required n-point 
correlation functions coming from this new representation are derived in appendices B 
and C. Finally, in section 4 we conclude. 



2. The contact process 

2.1. Field-theoretical description 

The master equation which describes the critical bosonic contact process as defined in 
section 1 can be turned into a field-theory in a standard fashion through an operator 
formalism which uses a particle annihilation operator a{r,t) and its conjugate a''"(r,t), 
see for instance [201 I22j for detailed discussion of the technique. Since we shall be 
interested in the connected correlator, we consider the shifted field and furthermore 
introduce the shifted response field 

0(r,t) := a{r,t) - po 

0(r, t) := a{r, t) = a\r, t) - 1 (10) 

such that (0(r, t)) = (our notation implies a mapping between operators and quantum 
fields, using the known equivalence between the operator formalism and the path-integral 
formulation [221 1221) ■ shall see, these fields and will become the natural 

quasiprimary fields from the point of view of local scale-invariance. We remark that the 
response function is not affected by this shift, since 

dh{r',s) oh[r\s) 
Then the field-theory action reads, where is the reaction rate j25] 
5[0,0]= jdRjdu ^{2Mdu - - fi^\<l) + Po) 

= S'o[0,0] + 540,0] (12) 

To keep expressions shorter, we have supressed the arguments of 0(-R, u) and 0(-R, u) 
under the integrals and we shall also do so often in what follows, if no ambiguity arises. 
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The diffusion constant D is related to the 'mass' A4 through D = {2Ai)~^. We have 
decomposed the action as follows: 

So[^,^]:= jdR jdu ^{2Mdu - (13) 

describes the deterministic,^ noiseless part whereas the noise is described by 

540, 0] := -fi [dR [ du \^\<f) + po)l . (14) 



quite analogously to what happens in the kinetics of simple magnets, see for details. 

In principle, an initial correlator G{r; 0, 0) could be assumed and will lead to a 
further contribution ^ini to the action. For critical systems, one usually employs a term 
of the form 5ini,st = JdR {(j){R, 0) - {(f){R, 0))f, see e.g. [HI E] but this would have 
for us the disadvantage that it explicitly breaks Galilei-invariance. We shall rather make 
use of the Galilei-invariance of the noiseless action 5*0 [0,0] and use as an initial term 

= -^j dRdR' 4>{R, {))G{R - R'- 0, O)0(i?', 0). (15) 

Because of the initial condition G(i?; 0,0) = 0, however, S'ini[0] = and we shall not 
need to consider it any further. 

From the action ()12|l . n-point functions can then be computed as usual 

(0i(r-i, ti) ... 0„(r„, tn)) = j VcpVcf) 0i(ri, ti) . . . 0„(r„, t„) exp (^-S[(f), 0] j (16) 

which through the decomposition ()12|) can be written as an average of the noiseless 
theory 

(0i(ri, ti) . . . 0„(r„, tn)) = <^0i(ri, ti) . . . 0„(r„, tn) exp 0]) (17) 

where (. . .)o denotes the expectation value with respect to the noiseless theory. 

2.2. Symmetries of the noiseless theory 

In what follows, we shall need some symmetry properties of the noiseless part described 
by the action So [0,0] which we now briefly recall. The noiseless equation of motion 
for the field is a free diffusion-equation 2Aidt4>{x,t) = V^0(a;,t). Its dynamical 
symmetry group is the well-known Schrodinger-group Sch{d) [221 HZI which acts on 
space-time coordinates {r,t) as {r,t) ^— {r',t') = g{r,t) where 

+ .S-,3,^l (18) 

^ This terminology is used since the equation of motion of (p foUowing from 6*0 is a partial differential 
equation and not a stochastic Langevin equation. 
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field 


scaling dimension 


mass 






X 


M 







X 


-M 






X2 


-2M 


T 


:=0V 


xr 


-M 


S 


:= 


Xy, 


-2M 


r 


= 0302 


Xy 


-M 



Table 3. Scaling dimensions and masses of some composite fields. 



and where 7^ is a rotation matrix. Solutions of the free diffusion equation are carried 
to other solutions of the same equation and transforms as 

0(r,t) (r,0)(r,t) = f,[g-\r,t)]<P[g-\r,t)] (19) 

where the companion function fg is known explicitely and contains the so-called 'mass' 
M = {2D)-^ EH 123. We hst the generators of the Lie algebra scf)i = Ue{Sch{l)) in 
one spatial dimension [221 

X-i = -dt 

1 X 

Xo = - tdt - -rdr - - 

Xi = - t^dt - trdr -xt- —r^ 

Y_i = -dr (20) 

Yi = -tdr-Mr 



2 



Mo = - M 

Fields transforming under Sch{d) are characterized by a scaling dimension and a mass. 
We list in table El some fields which we shall use below. We remark that for free fields 
one has 

X2 = 2x , xr = 2x + X , xs = 3x + X , xr = Sx + 2x (21) 

but these relations need no longer hold for interacting fields. On the other hand, from 
the Bargman superselection rules (see [30j and below) we expect that the masses of the 
composite fields as given in table El should remain valid for interacting fields as well. 

Throughout this paper, we shall make the important assumption that the fields 
and transform covariantly according to (fT^ under the Schrodinger group. By analogy 
with conformal invariance, such fields are called quasiprimary [8 . For quasiprimary 
fields the so-called Bargman superselection rules [SO] holds true which state that 

. . . )o = unless n = m (22) 



LSI: bosonic contact and pair-contact processes 



9 



We recall the proof of these in appendix B. Before we consider the consequences of ((221), 
we recall the well-known result on the form of noise-less n-point functions in ageing 
systems. 

Since in ageing phenomena, time-translation invariance is broken, we must consider 
the subalgebra age^ C sci)i obtained by leaving out the generator of time-translations 
X_i [31 . Then the n-point function of quasiprimary fields (pi, i = 1, . . .n has to satisfy 
the covariance conditons |29[ l8ij 



n 



(^i(ri,ti)...yp„(r„,t„))o = ; A; e {0, 1} (23) 



j=i 



f^ywj (^i(n,tO---^n(r„,t„))o =0 ; mG|-^,i| (24) 

where (fi stands either for a quasiprimary field (pi or a quasiprimary response field 
(pi. The (pi are characterized by their scaling dimension Xi and their mass Aii. The 
generators Xk are then the extension of (j^UI) to n-body operators and the superscript 
[i) refers to ipi. The n-point function is zero unless the sum of all masses vanishes 

n 

= (25) 

i=l 

which reproduces the Bargman superselection rule ()22|). It is well-known j^ lH] that the 
noiseless two-point function -Ro(^; t, s) = {fi{r, t)ip2{r, s))o is completely determined 
by the equations and up to a normalization constant. 

Mr, r'; t, s) = i?o(t, s) exp j S{Mi + M2) (26) 

where the autoresponse function is given by 

R,(t,s)=roit-s)-2^^^+^'^n-] (27) 



This reproduces the expected scaling form ((Zj) together with the scaling function fn^y) 
as given in table |21 if we identify 

X = Xi = X2, and x = a + 1 (28) 

For the critical bosonic contact process, we read off from tabled that a = f — 1- Hence 
one recovers x = |, as expected for a free-field theory. 



2. 3. Reduction formulcB 

We now show that the Bargman superselection rule ()22|1 implies a reduction of the n- 
point function of the full theory to certain correlators of the noiseless theory, which is 
described by Sq only. This can be done generalizing the arguments of [T^ . 



LSI: bosonic contact and pair-contact processes 10 

First, for the computation of the response function, we add the term 
JdR J du(j){R,u)h{R,u) to the action. As usual the response function is 



R{r,r';t,s) = ^0(r, t)0(r', s) 

{r,t)4>{r',s)exp(^-fi JdR j du^\R,u){^{R,u) + po)^'j 



{r,t)<p{r',s)) =Ro{r,r';t,s) (29) 
/ 

where we expanded the exponential and applied the Bargman superselection rule. 
Indeed, the two-time response is just given by the response of the (gaussian) noise- 
less theory. We have therefore reproduced the exact result of table El for the response 
function of the critical bosonic contact process. 

Second, we have for the correlator 

G{r,r',t,s)= ^0(r,t)0(r's)exp ^-/i JdR J du4^'^{R,u)^{R,u)^ 

X exp ^-/ipo jdRjdu 4>'^{R, u)^^ (30) 

Expanding both exponentials 

exp (^-fi JdR jdu 4>\R, u)4>{R, u)^ = ^ (/^'^ /^^ ^'>^^^' ^' 

exp (^-/ipo jdR jdu 4^\R, u)^ = ^"^f^ {/ J '^'^ 

and using the Bargman superselection rule ()22p. non- vanishing terms only arise if 
2n + 2m = + 2 or else 

n + 2m = 2 (31) 

This can only be satisfied for n = and m = 1 or for n = 2 and m = 0. The full noisy 
correlator hence is the sum of only two terms 

G(r, r'; t, s) = Gi(r, r'; t, s) + G^ir, r'; t, s) (32) 

where the first contribution involves a three-point function of the composite field 0^ of 
scaling dimension X2 (see table Ej) 

Gi{r, r'; t, s) = -/ipo JdRjdu (^0(r, t)0(r', s)^'^{R, u)^^ (33) 

whereas the second contribution comes from a four-point function and involves the 
composite field T (see table Ej) 

G2{r, r'; ^, s) = y J dRdR' J dudu' {(f){r, t)(j){r', s)T{R, u)T{R', m'))o (34) 
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We see that the connected correlator is determined by three- and four-point functions of 
the noiseless theory. We now use the symmetries of that noise-less theory to determine 
the two-, three,- and four-point functions as far as possible. 



2.4- Correlator with noise 

We consider Gi{r,r',t,s) first. The appropriate three-point function is given in 
appendix B, equation ()B24|) : 

(0(r, t)(P{r', s)0'(i?, u))o = it- (t - u)-^^'^'' (s - m)-^-^ 

X exp _ Mu„v,)e{t - u)e{s - u) (35) 

V 2 t — u 2 s — u 



with 



_u [{s - u){r - R) - {t - u){r' - R)f 



t (t — u)(s — u) 



2 



^^^u_ [(s-u){r-R)-{t-u){r'-R)]^ ^^^^ 

s (t — uy{s — u) 

and an undetermined scaling function \&3. The B-functions have been introduced 
by hand because of causahty but this could be justified through a more elaborate 
argument along the lines of [HI]. Introduced into ((221), ^^is gives the general form for 
the contribution Gi(r, r'; t, s). We concentrate here on the autocorrelator, i.e. r = r' 
and find, with y = t/s 

\y{y-e){i-er {y-eni-e)) ^'^^ 

where H is an undetermined scaling function. Very much in the same way, we find for 

G2{t, s) 

2 rl n 

2 Jo Jo 

X (, - - eri" Ij^R' exp (-^^ - 

X ^4 {usiR, e, R', e'), u^{R, e, r', e'), v3{r, e, r', e'),Vi{R, e, r', e')) (38) 

where \l/4 is another undetermined function and the functions U3,U4,V3, V4 can be worked 



out from the appropriate expressions ()B26|1 in the appendix B by the replacements 
rs — r2 — R, r4 — r2 ^ R', t2 ^ I, ti ^ y, ts ^ 9, t4^ ^ 6' (remember that ri = 

As we have a free-field theory for the critical bosonic contact process, we expect 
from table HI and eq. that x = x = d/2 and hence the following scaling dimensions 
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for the composite fields 

3 

X2 = d, xr = -d (39) 

Consequently, the autocorrelator takes the general form 

G{t, s) = s'-'l^gr{t/s) + s^''g,{t/s) (40) 

For d larger than the lower critical dimension (i* = 2, the second term merely furnishes a 
finite-time correction. On the other hand, for d < d^ = 2, it would be the dominant one 
and we can only achieve agreement with the known exact result if we assume \E'4 = 0. In 
what follows, we shall discard the scaling function g2 and shall concentrate on showing 
that our expressions for gi are compatible with the exact results given in table |21 

In order to do so, we choose the following special form for the function \E'3 

*3(wi,t^i)=Hf---) (41) 

where S remains an arbitrary function. Then we are back in the case already treated in 
Hg. We find 

G,{t,s) = -/iposi+i-^-^^^(i/- l)^"-i 

X ^'d^ [{y -e){l- 9)]i-h'-^, i^-^^~r) ^^^^ 



where the function 0i is defined by 

(Pi{w)= [dRexp(-^R'')E{R') (43) 



As in jTz] we choose 

Mw) = 0o,cw;-^-^ (44) 

This form for 0i(w) guarantees that the three-point response function 
(0(r, t)(p{r, s)0^(r', u))o is nonsingular for t = s. We have thus 

G{t,s) = G,{t,s) = s-'fG(-] (45) 



fciy) = - m0o,c / d9iy + l- 29)-l 
Jo 



with 



^ ^iy-ir^^'-{y + l)-^^') (46) 

and we have reproduced the corresponding entry in table |21 for the critical bosonic 
contact process."*" 

+ We remark that for 2 < d < 4, the same form of the autocorrelation function is also found in the 
critical voter- model jH2j . 
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3.1. Field-theoretical description and reduction formula 



For the pair-contact process we have two different cases, namely the case a < ac and 
the case at criticality a = ac- The following considerations apply to both cases and we 
shall for the moment leave the value of a arbitrary and only fix it at a later state. 

The action for the pair-contact process on the critical line is [23 eq. (30)] 



S[a, a] 



j dR j du [a{2Mdt - V^)a - aa^a^ - fxa^a^] 



(47) 



As before, see eq. (fTnjl . we switch to the quasiprimary fields 4>{r,t) = a{r,t) — po and 
4>{r,t) = a{r,t). Then the action becomes 



^[0,0] 



dR / du 



(l){2Mdt - V^)0 - acf^cj? - aplcj)'^- 



(4^ 



Also in this model, similarly to the treatment of section 2, a decomposition of the action 
into a first term with a non-trivial dynamic symmetry and a remaining noise term is 
sought such that the correlators and responses can be reexpressed in terms of certain 
n-point functions which only depend on 5*0. The first term reads 



5*0 [0,0] := 



dr / dt 



^{2Mdt - V2)0 - ai 



/,2 i2 



(49) 



and we derive its Schrodinger-invariance in appendix A. The remaining part is the noise- 
term which reads 



dR / du 



2 J,2 



2apo 



3 J,2 



2 J,3 



(50) 



Also in this case the Bargman superselection rule holds true. This means that we 
can proceed now in a very similar way as before.* First we have to check which n-point 
functions contribute to the response and correlation funtion. We rewrite exp(— 5'b[0, 0]) 
as a product of five exponentials and expand each factor. The indices of the sums are 
denoted by ki for the z-th term in (jSUj) . for instance for the first term 



exp 



dR / du 



2 J,2, 



R,u] 



oo ^ 

Y^ — i- IdR Idu 

fci=0 ^ 



2 J,2, 



R,u) 



(51) 



For the response function again only the first term of each sum contributes, that is 

Rir,r';t,s)=Roir,r';t,s) (52) 

* This argument works provided each term in Sb contains at least one response field (f> more than 
order-parameter fields (f>. 
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contr. 


IV I 


ko 






comp. field 


scaling dim. 


3-point /4-point 


Gi{t,s) 


1 















3-point 


G2it,s) 





2 








T 


xr 


4-point 


Gsit,s) 








2 





r 


xr 


4-point 


G^{t, s) 











1 


s 




3-point 


G^it, s) 





1 


1 





T and T 


xr,xr 


4-point 



Table 4. Contributions to the correlation function: The first column shows how we 
denote the contribution, the next four columns give the value of the corresponding 
indices. The sixth column lists the composite field(s) involved, the seventh column 
how we denote the scaling dimension of that field. The last column lists whether it is 
a three- or four-point function that contributes. 



is noise- independent. For the correlation function, we have the condition 2ki + 2/c2 + 
Sks + Sk^ + 3k^ = 2 + k2 + 2A;3 -|- k^ or simply 

2ki + k2 + k3 + 2^4 + 3A;5 = 2 (53) 

which implies immediately that 

k, = 0. (54) 

In table 0] we list the five differenent contributions to the correlation function. We 
denote also the form of the composite field, its scaling dimension and whether it is a 
three- or four-point function which contributes. A short inspection of the general form 
of the n-points function given in the appendix B shows that the contributions have the 
form (with y = t/s) 

Gi(t, s) = s---^^^+i+Vi(2/) , G,{t, s) = ^--+1+1/4(2/) (55) 
for the 3-point functions and 

G2it, s) = -^+'^+2/2(2/) , G^it, s) = s-™+'^+2/3(2/) 

G5(t,s)=S---5-T-|-'i^+'^+V5(2/) 

for the four-point functions. The scaling functions fi{y) involve an arbitrary functions 
\E'j which are not fixed by the symmetries (see appendix B for details). As we do not have 
a free-field theory in this case we can not make any asumptions about the value of the 
scaling dimensions of the composite fields. Therefore we do not know which terms will 
be the leading ones in the scaling regime. However, it turns out that the term Gi{t, s) 
alone can reproduce our result correctly. Thus we set the scaling functions /n=0 with 
n = 2, . . .5 analogously to the last section. We now concentrate on 

Gi{t, s) = apl dR du (^(f){r, t)(j){r, s)4>'^{R, u)^ ^ (56) 
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3.2. Symmetries of the noiseless theory 

As in the last chapter, we require for the calculation of the two- and three-point functions 
the symmetries of the following non-linear 'Schrodinger equation' obtained from 

2Mdt(j){x, t) = V^(j){x, t) + ^(0, 0) (57) 

with a nonlinear potential 

J^{(j),^) = -g(j)\x,t)^{x,t) (58) 

While for a constant g the symmetries of this equation are well-known, it was pointed out 
recently that g rather should be considered as a dimensionful quantity and hence should 
transform under local scale-transformations as well ^Tj. This requires an extension of 
the generators used so far and we shall give this in appendix A. The computation of 
the n-point functions covariant with respect to these new generators is given in the 
appendices B and C. In doing so, we have for technical simplicity assumed that to each 
field (fi there is one associated coupling constant gi and only at the end, we let 

g^ = . . . = g^=: g (59) 
Therefore, from eq. ()52|1 we find for the response function (see ()(7l2|) ) 



Ro{r,r';t,s) = {t - s)-^ 



M(r - r'Y\ ~ f t t- s q 
X exp - ■ -TAT. TT^ (60) 



t-s ) \s g^/y ' {t - s)y ^ 

with an undetermined scaling function \&2- This form is cleary consistent with our 
results in table |21 if we identify 

d ~ 

X := Xi = X2 = a + 1 = - , ^2 = const. (61) 

This holds true for both a < ac and a = ac- In distinction with the bosonic contact 
process, the symmetries of the noiseless part 5*0 do not fix the response function 
completely but leave a certain degree of flexibility in form of the scaling function \E'2- 

For the calculation of the correlator we need from eq. ()56|1 the following three-point 
function 

(r, t)0(r', s)4)\R, u)^ = {t - sy~^^^{t - u)-^^^{s - u)-^^^ 

X exp (-^<^ - MiJ^^ .A., .A) (62) 
\ 2 t — u 2 s — u I 



with 



_ u [{s-u){r-R)-{t-u){r' -R)f 
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_u [{s - u){r - R) - [t - u){r' - R)f 
^ s (t — uY{s — u) 



1 a^/y 
$2 [t - uy 



/?2 = - 
S2 



a 



u 



/?3 



S2 



1 1 

t — u u 



We choose the following realisation for 



{-- 









ia-b) 



16 

(64) 

(65) 
(66) 

(67) 



where the scaling function S was already encountered in eq. (j^T|) for the bosonic contact 
process. We now have to distinguish the two different cases a < ac and a = ac- For 
the first case a < ac, we have a — 6 = so that the last factor in ^7\i disappears and 
we simply return to the expressions already found for the bosonic contact process, in 
agreement with the known exact results. However, at the multicritical point a = ac 
we have a — b and the last factor becomes important. We point out that only the 
presence or absence of this factor distinguishes the cases a < ac and a = ac- 

If we substitute the values for Pi,P2 and Ps, ^3 becomes 



^3(Ml,t^l,/?l,/52,/?3) = S 



(-- 









0{y 



1 (a-b) 



{y-e){i-e) 



{61 



This factor does not involve R so that we obtain in a similar way as before 

."1 

-bf„, 1 \ {b—a)—a—l 



Gi(t,s) 



'y + l 



X 01 

where we have identified 

X2 = 2(& — a) + d 



29 



y-i 



%-i) 



{y-e){i-e) 



(69) 



(70) 



Gi(t,s) reduces to the expression Q if we choose the same expression for (pi{w) as 
before. We have thus reproduced all scaling functions correctly. 



4. Conclusions 

The objective of our investigation has been to test further the recent proposal of using 
the non-trivial dynamical symmetries of a part of the Langevin equation in order to 
derive properties of the full stochastic non-equilibrium model. To this end, we have 
compared the known exact results for the two-time autoresponse and autocorrelation 
functions in two specific models, see table with the expressions derived from the 
standard field-theoretical actions which are habitually used to describe these systems. 
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This is achieved through a decomposition of the action into two parts S = Sq + Sf, such 
that (i) 5*0 is Schrodinger-invariant and the Bargman superselection rules hold for the 
averages calculated with So only and (ii) the remaining terms contained in Sb are such 
that a perturbative expansion terminates at a finite order, again due to the Bargman 
superselection rules. The two models we considered, namely the bosonic variants of the 
critical contact and pair-contact processes, satisfy these requirements and are clearly 
in agreement with the predictions of local scale- invariance (LSI). In particular, our 
identification eq. of the correct quasi-primary order-parameter and reponse fields 
is likely to be useful in more general systems. 

Specifically, we have seen the following. 

(i) In the bosonic contact process, the symmetries of the noiseless part 5*0 of the action 
is described in terms of the representation of the Schrodinger-group relevant for the 
free diffusion equation. 

In consequence, the form of the two-time response function is completely fixed by 
LSI and in agreement with the known exact result. The connected autocorrelator is 
exactly reducible to certain noiseless three- and four-point functions. Schrodinger- 
invariance alone cannot determine these but the remaining free scaling functions 
can be chosen such that the known exact results can be reproduced. 

(ii) For the bosonic pair-contact process, the symmetries of the partial action 5*0 are 
described in terms of a new representation pertinent to a non-linear Schrodinger 
equation. This new representation, which we have explicitly constructed, involves a 
dimensionful coupling constant g. Therefore even the response function is no longer 
fully determined. As for the autocorrelation function, which again can be exactly 
reduced to certain three- and four-point functions calculable from the action ^o, 
the remaining free scaling functions can be chosen as to fully reproduce the known 
exact results. 

The consistency of the predictions of LSI with the exact results of these models furnishes 
further evidence in favour of an extension of the well-known dynamical scaling towards 
a (hidden) local scale-invariance which influences the long-time behaviour of slowly 
relaxing systems. An essential ingredient were the Bargman superselection rules which 
at present can only be derived for a dynamical exponent z = 2. An extension of our 
method to models with z 2 would first of all require a way to generalize the Bargman 
superselection rules. We hope to return to this open problem elsewhere. 
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Appendix A. Representations of age^ and sc()^ for semi-linear Schrodinger 
equations 

We discuss the Schrodinger-invariance of semi-linear Schrodinger equations of the form 
fl57|) and especially with non-linearities of the form (j^Hjl . With respect to the well- 
known Schrodinger-invariance of the linear Schrodinger equation, the main difference 
comes from the presence of a dimensionful coupling constant g of the non-linear term. 

It is enough to consider explicitly the one-dimensional case which simplifies the 
notation. In one spatial dimension, the Schrodinger algebra 5ci)i is spanned by the 
following generators 

5C[)i = (X_i,Xo,Xi,r_i/2,>^l/2,Mo> (Al) 

while its subalgebra agc^ is spannned by 

agei = (Xq, X,, Y^,/2, Fi/a, Mq) (A2) 

These generators for g = are listed explicitly in eq. ()2Up and the non-vanishing 
commutators can be written compactly 

[Xn, Xn'] = (n — n')Xn+n' 

[Xn, Ym] = in/2- m)Yn+m 

[Fi,F.]=Mo (A3) 

where n,n' E {±1,0} and m E {±|} (see for generalizations to d > 1). 

Following the procedure given in ^7j, we now construct new representations of agc^ 
and of sci)i which takes into account a dimensionful coupling g with scaling dimension 
y as follows. 

(i) The generator of space-translations reads simply 

Y_i = -dr. (A4) 

(ii) The generator of scaling transformations is assumed to take the form 

1 X 

Xq = -tdt - -rdr - ygdg - - (A5) 

where y is the scaling dimension of the coupling g. 

(iii) For sz\)i we also keep the usual generator of time-translations 

X_i = -dt. (A6) 

(iv) The remaining generators we write in the most general form adding a possible 
^f-dependence through yet unknown functions L, Q, P. 

Mo= - M- L{t,r,g)dg 

Yi = - tdr- Mr -Q{t,r,g)dg (A7) 



Xi = - t'^dt - trdr - - xt - P{t, r, g)d. 



9 
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The representation given by eqs. ()A4IA5IA6IA7|) must satisfy the commutation relations 
()A3|) for agti or sci)^. From these conditions the undetermined functions L,Q and P 
are derived. A straigthforward but shghtly longish calculation along the lines of [T7] 
shows that for ogc^, one has 

L = , g = , P = po{M)ty+^m{t/g) (A8) 

Here, m{v) is an arbitrary differentiable function and pQ{Ai) a 7W-dependent constant. 
We shall use the shorthand v = t^ / g in what follows. 

In consequence, for oge^ only the generator Xi is modified with respect to the 
representation eq. ()2()j) and this is described in by the function m{v) and the constant 
Po{M). 

On the other hand, for scf)^ the additional condition = 2Xq leads to 

Po = 2y, m{y) = v~^. 

Hence, the new representations are still given by eq. ()2()j) with the only exception 
of Xi which reads 

agei : Xi = -t^dt - trdr - po{M)ty+^m {t^ g) dg xt 

sc[)i : Xi = -t^dt - trdr - ^ytgdg xt (A9) 

We require in addition the invariance of linear Schrodinger equation {2Aidt — d^)(f) = 
with respect to this new representation. In terms of the Schrodinger operator S this 
means 

[S, X] = XS ; where S := 2MoX_i - Y^^^^ (AlO) 

and X is one of the generators of age^ eq. ()A2|) or of sci)i eq. ()A1|) . Obviously, A = 
if A" G {X_i,Y±i/2, Mq) and A = — 1 if A" = Xq. Finally, for Xi we have from the 
definition of the Schrodinger operator S 

S, Xi] = - AMoXo + {Y1/2Y-1/2 + y-1/211/2) 

= - 2tS + M {1 - 2x - Aygdg) (All) 

where in the second line the explicit forms eqs. ()A4|A5jA6|A7j) were used. This also 
holds for all those representations of agc^ for which there exists an operator X_i ^ aQti 
such that = 2Xo and we shall restrict our attention to those in what follows. 

On the other hand, the direct calculation of the same commutator with the explicit form 
gives for age^ 



S,Xi 



-2tS + M{l~2x)- 2Mpo{M)ty [{y + l)m{v) + yvm'{v))] dg (A12) 



Besides A = —2t, the consistency between these two implies for m{v) the equation 

v({y + l)m{v) + yv^^) = ^ (A13) 
V dv J po 
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with the general solution 



m(v) 



PO V Po 



(A14) 



where mo = mo{Ai) is an arbitrary constant. The larger algebra 5ci)i is recovered from 
this if we set po = 2y and mo = 0. Hence the final form for the generator Xi in the 
special class of representations of the algebra aQCi defined above is 



Xi = -t^dt - trdt - 2ytgdg - mog 



xt 



(A15) 



Summarizing, this class of representations of agc^ we constructed is characterized by 
the triplet (x, A^,mo), whereas for sci)i, the same triplet is {x,A4,0). 

Finally, to make Xi a dynamical symmetry on the solutions $ = $g(t, r) of the 
Schrodinger equation S^g = we must impose the auxiliary condition (1 — 2x — 
Aygdg)(^g = which leads to 



'^9it^r) = g 



il-2x)/{4y) 



$(t,r) 



(A16) 



In particular, we see that if x = 1/2, we have a representation of agCi without any 
further auxiliary condition. 

We now look for those semi-linear Schrodinger equations of the form 
5$ = F(t, r, (jf, $, $*) for which the representations of oge^ or sc[)i as given by 
eqs. flA4jA5|A6|A7|l and with Xi as in ()A15|1 act as a dynamical symmetry. The non- 
linear potential F is known to satisfy certain differential equations which can be found 
using standard methods, see [23,1111 eq. (2.8)]. In our case these equations read 



Y_i 
Mo 
Yi 

2 



dtF = 
drF = 

($9$ - - 1)F = 

[tdrF - Mr{^d^ - ^*d^* - 1)] F = 

tdt + ^rdr + ygdg + 1 - |(<l'9$ + <l>*d^* - 1) 
t^dt + trdr + 2t{ygdg + 1) + niog^+^^ydg 



(A17) 
(A18) 
(A19) 
(A20) 

(A21) 



F = (A22) 



We first solve these for sc()i. From the conditions eqs. ()A17IA18IA19IA20IA2T|) we 

easily find 



F = $(<l>$*)^/^/(^^($$*)^ 



(A23) 



where / is an arbitrary differentiable function. Two comments are in order: 
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(i) For a dimensionless coupling g, that is y = 0, we have x = 1/2. Then the scahng 
function reduces to a (^-dependent constant and we recover the standard form for 
the non-hnear potential F as quoted ubiquitously in the mathematical literature, 
see e.g. jSl]. 

(ii) Taking into account the generator Xi from eq. ()A22|) as well does not change 
the result. Hence in this case translation-, dilatation- and Galilei-invariance are 
indeed sufficient for the special Schrodinger-invariance generated by Xi, see also 
|31j . We point out that traditionnally an analogous assertion holds for conformal 
field-theory, see e.g. |1 , but counterexamples are known where in local theories 
scale- and translation-invariance are not sufficient for conformal invariance 



Second, we now consider the representation of age^ where Xi is given by ()A15|1 . 
We have the conditions eqs. (IA18IA19IA20IA21IA22ll . We write F = ^T{uj,t,g) with 
LJ = $$* and the remaining equations coming from Xq and Xi are 

{tdt + ygdg - xudu + 1) = 

[eOt + mog^+^/ydg) T = (A24) 
with the final result 



F = $($$*)'/^/ ($$ 



g-yy _ 

yt. 



(A25) 



and where / is the same scaling function as encountered before for scf)]^. Finally, the 
result for the general representations of age^ which depend on an arbitrary function 
m{v) are not particularly inspiring and will not be detailed here. We observe 

(i) For mo = 0, this result is identical to the one found for sci)^. 

(ii) Even for mo 7^ 0, the form of the non-linear potential reduces in the long-time limit 
t ^ 00 to the one found in eq. ()A23j) for the larger algebra sci)i. 

We can summarize the main results of this appendix as follows. 



Proposition. Consider the following generators 

Mo = - M , F_i/2 = -dr , Yi/2 = -tdr " Mr , X_i = -dt 
Xo= - tdt - \rdr - ygdg - ^ (A26) 

Xi= - t^dt - trdt - 2ytgdg - m,g^+^'ydg - ^ - xt 

where x,Ai,mQ are parameters. Define the Schrddinger operator S := 2MoX_i — Y^-^^^^. 
Then: 

(i) the generators (Xo,i, Ki-1/2, Mq) form a representation of the Lie algebra agej. If 
furthermore rriQ = 0, then (Xo,±i, l±i/2. Mo) is a representation of the Lie algebra sci)i. 

(ii) These representations are dynamical symmetries of the Schrddinger equation = 
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0, under the auxiliary condition (1 — 2x — Aygdg)^ = 0. 

(Hi) For the Schrddinger- algebra scf)^ and also in the asymptotic limit t ^ oo for 
the ageing algebra age^, the semi-linear Schrddinger equation invariant under these 
representations has the form 

5$ = $ f (^gx (A27) 

where f is an arbitrary differentiate function. 

This general form includes our potential ()58|) since the scaling dimension y is a 
remaining free parameter in our considerations 



Appendix B. The n-point function 

In this appendix we use the generators ()A26|) from appendix A to find the most general 
form of the n-point functions compatible with the symmetries for > 3. We shall do 
this for the case mg = only, as this will be enough to reproduce the exact results of 
table El The case n = 2 needs a special treatment and is presented in appendix C. 

We restrict ourselves to the case (i = 1 for simplicity, but the generalisation to 
arbitrary dimension will be obvious. First we introduce some notation. We fix an 
arbitrary index k and define the shifted coordinates 

f-b := Vb - Vk, ib := tb - 4 for j ^ k and 4 := tk (Bl) 

In the sequel, we will adopt the following convention: The index a always runs from 1 
to n, the index b runs from 1 to n but skips k. The prime on a sum means that the 
index k is left out, viz. 

n ' n 

E^-=E^' (B2) 

i = l i=l 

We denote the ?7,-point function by 

FiWa}, {ta}, {ga}) ■= {fi{ri,ti) . . . (/?„(r„, t„)) (B3) 

where we assume one coupling constant for each field. This quantity has to satisfy the 
following four linear partial differential equations 

E^?Mi^(K},{U,{U) = o ; ke{0,l} (B4) 
i=l / 

{Y.J^^^^(M,{ta},{ta}) =0 ; (B5) 
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To solve these equations, we use the method of characteristics |27j. We solve ()B4|) first 
for spatial translation invariance with the result 

F{{ra}, {ta}, {ga}) = F{{h}, {ta}, {9a}) (B6) 

with a new function F with 3n — 1 arguments. In order to solve for Xq we set 
1 " 

x = -J2^^ (B7) 



2 

i=l 



and make the ansatz 

Hih}, {ta}, {ga}) = lliU - tj)-''^^Gi{h}, {ta}, {^7a}) (B8) 



i<j 



where the parameters pij and the function G remain to be determined. We also change 
to the new independent temporal variables ia- Then one finds after a short calculation 

' n 1 " ' ^ \ 

+ 2E +J2y^9^^g, ^9a}) = (B9) 

V i=l i=l 1=1 / 

together with the condition 

J2P^r (BIO) 



X 

i<j 



Before proceeding to solve this equation, we turn to the generators li/2- We find for 
the function G{{fb}, {ia}, {ga}) 

' n ' n ' \\ 

iidf, + -^^^^ + 5^ -^O {4}, {go}) = 0. (Bll) 

vi=l i=l \i=l / / 

Since G{{rb}, {ia}, {ga}) does not depend on r^, we recover the Bargman superselection 
rule 

n 

YMi = (B12) 

i=l 

as expected. For G{{rb}, {ia}, {ga}) we make another ansatz: 

G{{h}, {ia}, {ga}) = exp ^^^^^ ^^^3) 



.=1 2 U 



where the function H{{rh}, {ta}, {ga}) remains to be determined. With ()B12|) and 
flB13|) . equation ()B11|) reduces to 

Y UdfA H{{h}, {ia}, {ga}) = 0. (B14) 
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We retake ()B9|) and introduce the ansatz ()B13|) . This yields 

(n ^ n ' n \ 

+ 2E ^^^^^ + Y.y^A H{{h}, {tj, {g,}) = (B15) 
i=l i=l 1=1 J 



The last equation for H{{rb}, {ta}, {ga}) we obtain from Xi. Using the ansatz (jH 
yields an equation for G 

1=1 i=l i=l i=l \i=l i=l J 

(n ' 1"' \\ 

Y + ""^^^^ + Y y^9^^9. Giin}, {ia}, {g,}) = 0. (B16) 
1=1 i=l i=l / / 

Together with the condition 

n 

+ = (B17) 

i<j i=l 

which is satisfied if we choose the parameters pij such that 

Xl = P12 + Pl3 + Pl4 + Pl5 + • • • + Pin 
X2 = P12 + P23 + P24 + P25 + • • • + P2n 

3^3 = Pl3 + P23 + P34 + P35 + • • • + P3n (B18) 



Xn — pin + P2n + P3n + P4?i + • ■ ■ + Pn-1 n 

Here a few remarks are in order. The above system is compatible with ()B10|) . as can be 
see be adding all equations. Also, this system is always solvable for n > 3, as for n > 4, 
it is underdetermined and for n = 3 the correponding determinant does not vanish jj. 
Lastly, we often have the case xi = X2 =: x and X3 = . . . = Xn ='■ x. In this case, we 
can set 

P12 = X - ^x ; p2i = \x, pii = \x for i = 3, . . . , n (B19) 

and Pij = for all the remaining pij. We still have to rewrite equation ()B16|1 in terms 
of the variables {vh} and {ta}- Here we take equations ()B9|) and ()B11|) and the ansatz 
()B13|) into account and get for H{{fb}, {ia}, {ga}) 

t\di^ - ildt, + Y i^^id^^ + 2^ y^^^Sid,^ ^(i^'-}' ^9a}) = 0. (B20) 

\i=l i=l i=l J 

We thus have to solve the homogenous equations ()B14|) . ()B15|) and ()B20|) . This will 
eliminate three more variables and yields 

F{{ra}. {ta}. {ga}) = ' t,)'^'^ exp (^-1^ ^^^^^^J^^^ 

tt This system is not solvable for n = 2 when xi ^ X2. This case is considered in appendix C. 
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x^„(K},K},{/3J) (B21) 

with an arbitrary function which depends on 3n — 4 variables. Here the index c runs 
from 1 to n but skips k and another arbitrarily fixed index r ^ k, and the expressions 
Uc, Vc and Pa are given by 

^ _ tk {{vc - rk){tr - tk) - (vr - rk){tc - tk) f 

(tc — tk)(tr — tk^tc 

tk {{rc - rk){tr - tk) - {Vr - rk){tc - tk)f 
= g^iVv.) ( , (3, = g,Wy^^ ( th^^LZlf.) (B22) 

V [tr — tk )tk J \ [tr — tk )tr J 

We remind the reader of our convention that the index c runs from 1 to n skipping r 
and k and that the index h runs from 1 to n skipping only k. 

In higher dimensions rotational invariance has to be satisfied as well and then the 
generalization to arbitrary d is straightforward. 

If we consider instead the algebra oge^ with dimensionless couplings gi we merely 
have to make the replacement 

^„ ({«c}, {Vc}, {(3a}) ({«c}, {Vc}) (B23) 

where \Ef„ is also an arbitrary function such that only the dependence on {(3a} drops 
out. 

Finally, we explicitely list the three- and four-point functions in the form in which 
they are needed in the main text. The three-point function with fixed indices r = 2 and 
/c = 3 and the special situation ()B19|) assumed reads 

F{{ra}, {ta}, {ga}) = (^1 - t^)~^'^-"^''\tr - t,)-"^'\t2 - ts)"^' (B24) 



with 



tl {tl - t3){t2 - t^) 



2 



ts [{ri - r3)(t2 - ts) - (r2 - r,)it^ - ts)f 



The four point function with r = 1 and = 2 in the special situation ()B19|1 reads 

F{{ra}, {ta}, {9a}) = {tl " t,)-(^-'\t^ - t,)-'^%t, - t,)-"^%t2 - ts)"^' 



X (t2-t4)-^"exp f -^E-^ /''; I ^n(K},M,R}) 
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with 

U3 



M4 



V4, 



to [frq 


_ r')){t^ — 




— r-yMt-^ 


-tA]? 






'-2 J l,'-! 






t2 [(r-4 


-r-2)(ti- 


t2) - (ri 


-r2)(t4 








-t2)(ti- 


t2)2 




t2 [(r-3 


-r-2)(ti- 


t2) - (r-i 


-r2)(t3 


-t2)f 


tl 




-t2)(t3- 


t2)2 




t2 [{U 


-r-2)(ti- 


t2) - (r-i 


-r2)(t4 


-t2)f 


tl 




-t2)(t4- 


t2? 





^3 = z u TUl TT^ (^26) 



My: 



'2 



Appendix C. The two-point function 

In this appendix we calculate the two-point function, which was not included in the 
treatment of appendix B. Again, we only treat the case mo = 0. Apart from the 
generator Ai, the calculations are similar to those done in appendix B, so we only give 
the essential steps. First we define 

T ■=ti-t2 , r := ri - r2 (CI) 

and then we proceed as follows. We solve for Mq, F_i/2, ^1/2, -^0 in exactly the same 
way as before with the result 

F{ri, r2, tl, t2, 9i, 92) = {^liri, ti, 5-1)^52(^2, ^2, 92))o = ^~''G'(r, r, t2, 91, 92) (C2) 
where x = ^ {xi + X2) and G{r, r, t2, 91, 92) satisfies the equations 

(^dr + t2dt., + ]^rdr + yi9idg^ + y292dg^ G{r, r, t2, 91, 92) = (C3) 

iTdr + rMi)Gir,T,t2,9u92) =0 (C4) 

and the Bargman superselection rule 

Mi + M2 = (C5) 

holds true. Now ()C3|) is solved by 

G{r, T, t2, 9i, 92) = G{ui, U2, Vi,V2) (C6) 

where we have defined 

^2 „2 „i/^'2 

' ' 9i 92 /r^^\ 

Ui := — , U2 := — , Vi := , V2 := (C7) 

T t2 T T 
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and rewriting ()C4|) in terms of the new variables yields 

uidu^ + M2<9„2 + ^uiMi^ G{ui, U2, Vl, V2) = (C8) 
which is solved by 

G{ui,U2,vi,V2) = exp I -]-uiMi] H{w,Vi,V2) , w:=— (C9) 

The function H{w, Vi,V2) is found through the generator Xi. Using again the invariance 
under Y1/2 and Xq, we readily obtain in terms of f 1, f2 and w 

[w + l)dw + vidy^ - V2d^2 + ^(xi - X2)^ H{w, vi,V2) = 0. (CIO) 
The most general solution of this equation is 

H = iw + l)-^(-i--^)^2 (^^^^^' ^1^2) (Cll) 

where the function \I^2 remains arbitrary. Substituting back the values for Vi,V2 and w 
our final result is 

'+ \ -3(^1-^2) 
1 \ 



^('^i, h, r2, t2) = Sm^+M2,o {h - t 



2 



l(xi+X2) 



ti-t2 J \\hj 9i iti-h, 



For applications to semi-linear equations, one now sets g '■= gi = 92 with a scaling 
dimension y ■= yi = y2- In the limit y ^ 0, the function ^2 reduces to a (yf-dependent 
normalization constant and we recover the standard result 



In many applications, one expects the scaling functions to be universal, up to 
normalization. On the other hand, the coupling g should be a non-universal quantity 
so that a universal scaling function cannot contain g in its arguments. This leads to 
^2 = ^2((ii/^2)^) and we point out that such a scaling form would be compatible (one 
still has z = 2, however) with what is found from the field-theoretical renormalization 
group and numerical simulations in non-equilibrium critical dynamics El- An 
extension to different values of z would as a first step require the generalization of the 
Bargman superselection rules. We hope to come back elsewhere to this open problem. 
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